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Abstract 


We investigate the helicity dependent twist-two and twist-three generalized 
parton distributions in light-front Hamiltonian QCD for a massive dressed 
quark target. Working in the kinematical region £ < x < 1, we obtain the 
splitting functions for the evolutions of twist-two quark and gluon distribu¬ 
tions in a straightforward way. For the twist-three distribution, we find that 
all contributions are proportional to the quark mass and thus the twist-three 
distribution is directly related to the chiral symmetry breaking dynamics in 
light-front QCD. We also show that the off-forward Wandzura-Wilczek rela¬ 
tion is violated in perturbative QCD for a massive dressed quark.We calculate 
the quark mass correction to the WW relation in the off- forward case and 
show that it is related to h\ (x) in the forward limit. We extract the ’genuine 
twist three’ part of the matrix element in the forward case and verify the 
Burkhardt-Cottingham and Efremov-Leader-Teryaev sum rules. 

Keywords: Generalized parton distributions, Light-front Hamiltonian, Perturbation theory 
PACS : ll.10.Ef, 12.38.Bx, 13.60.Hb, 11.15.Bt 


I. INTRODUCTION 

The generalized parton distributions (GPDs) [1] are being studied intensively in recent 
years. GPDs are hybrids of the usual parton distributions measured in inclusive processes 
like deep inelastic scattering (DIS) and form factors measured in elastic exclusive processes. 
In general, they can be expressed as off-forward matrix elements of light cone bilocal oper¬ 
ators. In the forward limit, GPDs reduce to normal parton distributions, which can then 
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be expressed as forward matrix elements of light-cone bilocal operators. The moments of 
GPDs over the parton momentum fraction x give form factors, which are given in terms of 
off-forward matrix elements of local operators. Thus, GPDs have a much richer structure 
and they connect various processes, both inclusive and exclusive. They provide new and 
important informations about the structure of the hadron. They can be probed in deeply 
virtual Compton scattering (DVCS) and hard exclusive production of vector mesons (for 
recent reviews on GPDs and hard exclusive reactions, see [2-4] and references therein). 

The GPDs have been investigated recently in light-cone formalism by several authors 
and an overlap representation of the plus component in terms of light-cone wave func¬ 
tions has been given [5,6]. GPDs have also been constructed using light-cone model wave 
functions [7]. The transverse and the minus components are somewhat complicated since 
they involve the constraint held ?/?_. They are usually called ’bad’ components, since the 
operators in these cases involve interaction terms and are higher twist objects. In other 
words, they involve direct quark-gluon dynamics. To interpret the experimental results for 
DVCS on a proton target in the presently accessible Q 2 range [8-11], it is of primordial 
importance to understand the effect of higher twist components. The perpendicular, or the 
twist-three component of the operator has been investigated using the Wandzura-Wilczek 
approximation [12], where the explicit interaction dependent parts of the operator as well 
as the quark mass terms were neglected. In this case, the twist-three matrix element can 
be expressed solely in terms of twist-two GPDs. In the forward limit, these relations reduce 
to the Wandzura-Wilczek relation for the transversely polarized structure function gr [13]. 
There is no theoretical justification for this approximation. In the forward case, recent ex¬ 
perimental results for the transverse polarized structure function indicate that the deviation 
from the Wandzura-Wilczek approximated form is substantial in some kinematic range for 
a nucleon target [14]. Therefore, it is interesting to make a full calculation of the off-forward 
twist-three matrix element, within the context of perturbative QCD, taking into account 
the explicit interaction dependence of the operator, the mass as well as the intrinsic trans¬ 
verse momentum of the partons. A convenient tool is based on the light-front Hamiltonian 
description of composite systems utilizing many-body wave functions. Instead of a hadron 
target, here we consider a simpler target like a quark dressed with a gluon and calculate 
the off-forward matrix elements within the context of perturbative QCD. The two-particle 
wave function is given in terms of the light-front QCD Hamiltonian [15]. This approach has 
been used extensively in the recent years to calculate polarized and unpolarized distribution 
functions in DIS, twist-two [16,17], twist-three [18,19] and twist four [20], as well as the 
transversity distribution [21], and the transverse momentum dependent distributions [22], 
Recently we have also extended it to calculate the off-forward matrix elements of light-front 
bilocal vector operators [23]. We verified the hclicity sum rule in perturbation theory and 
showed the effect of quark mass in the twist-three matrix element, which is absent in the 
forward limit. 

In this work, we investigate the hclicity dependent generalized distributions, which are 
the off-forward matrix elements of light-front axial vector operators, in light-front Hamil¬ 
tonian perturbation theory for a dressed quark target. We restrict ourselves to the region 
f < x < 1, of the generalized distributions, where f is the skewedness. In this region, the 
contributions come from the overlaps of two-body wave functions upto 0(a s ). We obtain 
the splitting functions corresponding to the evolution of the twist-two helicity dependent 


2 



distributions in a straightforward way. In the twist-three distribution, we show the contri¬ 
butions from the quark-gluon interaction dependent part, the quark mass and the intrinsic 
transverse momentum dependent parts of the operator explicitly. We find that all the three 
contributions are proportional to quark mass. We find that the Wandzura-Wilczek (WW) 
relation is not satisfied in perturbative QCD in the off-forward case for a massive quark, 
analogously as its forward counterpart. Our results also show that the twist-three distri¬ 
bution is directly related to the dynamical effect of chiral symmetry breaking in light-front 
QCD. We calculate the mass correction to the WW relation in the off-forward case. This 
contribution is related to h\{x) in the forward limit. Using the quark mass correction to the 
WW relation, we also obtain the ’genuine twist three’ part of the matrix element. We show 
that the first and second moments of this are zero which give the Burkhardt-Cottingham 
(BC) [24] and Efremov-Leader-Teryaev (ELT) [25] sum rules respectively. 

The plan of the paper is as follows. In section II, we investigate the twist two hclicity 
dependent distributions, involving both the quark and the gluon operators, for a dressed 
quark state. In section III, we calculate the off-forward matrix element of the transverse 
component of the axial vector current. We investigate the WW relation in the off-forward 
case and the quark mass correction to it in section IV. Summary and discussions are given 
in section V. The light front spinors for longitudinally and transversely polarized quarks 
are given in Appendix A. An outline of the derivation of the quark mass term in WW 
approximation is given in Appendix B. 


II. HELICITY DEPENDENT TWIST-TWO DISTRIBUTIONS 

A. Quark Distribution 


The twist-two hclicity-dependent generalized distribution is given by 

Uy = / A- e j*P + ‘-(P'A' | ^(_1_ )7 + 7 ^ ( 1_) | px). (2.1) 

Here, P, P' are the 4-momentums and A, \' are the hclicities of the initial and final states 
respectively. 

We work in the so called symmetric frame [5,6]. The momentum of the initial state is P Al 
and that of the final state is P' tl . The average momentum between the initial and final state 
is then P^ = IHL + P - . The momentum transfer is given by A M = P'^ L — P M , P[_ = — P± = 
skewedness £ = Without any loss of generality, we take £ > 0. We also get A” = 

The above matrix element is conventionally parametrized in terms of the hclicity de¬ 
pendent distributions, H(x,£,t) and E(x,£,t), where t is the invariant momentum transfer. 
The matrix element can also be expressed in terms of overlaps of light-front wave functions. 
The operator is given by, 

/ ^ |,PV by)7Vlt(£) = T/ (2.2) 

Here ip + = A + ^, = ( 77 ° 7 ± . The above expression is given in the light-front gauge, A + = 

0, where the path-ordered exponential between the fermion fields in the bilocal operator is 
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unity. For simplicity we suppress the flavor indices. In the two-component representation, 
we have the dynamical fermion field, 


Efp) = J2x\ 

A 


dp + d 2 p L 


[b(p, \)e~ iqz 


+ d t (p, X)e ipz ], 


and the dynamical gauge field, 


A\Z) = £ 

A 


dq + d 2 q ± 
2 ( 2 vr) 3 g+ 


[ 6 la(g,A)e-^ 


+ h.c. , 


(2.3) 


(2.4) 


Here, %a is the eigenstate of a z in the two-component spinor of 0 + . We have used the 
light-front 7 matrix representation: 


7 


0 


0 -i\ 
i 0 ) ’ 


7° 


0 i \ 
i 0 ) ’ 


7 


-id* 0 \ 

0 id* J ’ 


(2.5) 


with d 1 = o ' 2 and d 2 = —a 1 . e*(A) is the polarization vector of the transverse gauge field. 
The Fock space expansion of the operator is given by, 


o +5 = 253 


f dk ,+ d 2 k' L 


' 0 2(27t) 3 a Jit ' 2(27r) 3 \//E 
A(2xP + - U + - A’ + )6 t (A:, «)&(&', s') 

+ 5(2xP + + U + + A: + )cZ(A:, -s)d\k', -s') 
+ d(2PP+ + k + - k’ + )d(k, -s)b(k', s') 

+ 5(2xP + + k ,+ - k + )b\k, s)d\k', -s') 


xt 3Xs'- 


( 2 . 6 ) 


We have, k + > 0 ,k' + > 0, k + — k' + = p + —p ,+ = 2£p + . In the kinematical region, £ < £ < 1, 
only the first term in Eq. (2.6) contributes [ 6 ]. We restrict ourselves to this kinematical 
region. 

We take the state | P, a) of momentum P and helicity a to be a dressed quark consisting 
of bare states of a quark and a quark plus a gluon: 


P,a)= b\P,cr) | 0) 

dkfd 2 kj 


E 

(71, A 2 


dkjd 2 kj 


^2(2?r) 3 P+d 3 (P - ifej - fc 2 ) 


2(27r) 3 A:i j y 2(27r) 3 /c2 

<h(P,o I *1. <7i; *2, A. 2 )M(A: 1 ,cri)a t (*:2,A2) | 0). 


(2.7) 


Here ad and are bare gluon and quark creation operators respectively and 0i and 02 are 
the multiparton wave functions. They are the probability amplitudes to find one bare quark 
and one quark plus gluon inside the dressed quark state respectively. Up to one loop, if one 
considers all kinematical regions, there will be non-vanishing contributions from the overlap 
of 3-particle and one particle sectors of the state, this situation is similar to QED [5]. In the 
kinematical region we are considering, such kind of overlaps are absent and it is sufficient 
to consider dressing only by a single gluon. The state is normalized as, 
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( 2 . 8 ) 


(P‘, \' | P,X) = 2(2jt ) 3 P + S x ,x'S(P + - P ,+ )6 1 (P- L 

0i actually gives the normalization constant of the state [16]: 

i 12 1 a s f1 + x 2 Q 2 

0i =1-7—6/ / dx--<03—, 

2ir io 1 — x /// 

within order a s . Here e is a small cutoff on x. 

The matrix element becomes, 


P 


'-L, 


F + = 


y/l - - x) 


Y / ^Q-lWsuX 


S 1,52,A ' 


x -T 
!-<£ 


, + A± )^i fr3 ^2^ 2 A 


,X + ^ 

'1 + r 1 


(2.9) 


We have introduced Jacobi momenta Xj,^ -1 such that X)»Xj = 1 and 
boost invariant wave functions, 


( 2 . 10 ) 

0 and the 


0i = 0i, 


V’2(Zi,® L ) = V^ + 0(^ + ,^). 


( 2 . 11 ) 


The hrst term in Eq. (2.10) is the contribution from the single particle sector and the 
second term is the contribution of the two-particle sector of the state. Using light-front 
QCD Hamiltonian, the two-particle wave function is given in terms of 0i as : 


«U*.9 X ) = -hWUr- 


9 t 

- xl. 


(l 1 ) 2 \/(l-i) 72(2 tt)3 

(!-x) 


o ( r , ^ w ~± 

2--T-<7 


imo 


x 


L 1 — x 

WU001- 


x 


( 2 . 12 ) 


where g is the coupling constant, T a is the usual (| of Gell-Mann) color matrix and m is 
the bare mass of the quark. In the denominator of the above expression, we have neglected 
terms of order m 2 compared to (q L ) 2 ■ Using Eq. (2.1) and Eq. (2.12), we see that the 
linear mass terms which cause helicity flip, are suppressed in the matrix element. The terms 
quadratic in mass do not flip helicity, but they are suppressed too. We calculate the helicity 
non-flip part. 

Using Eq. (2.12) we get, 


F + 


i-i 2 < 5(1 


x) + 


2n 


Cflog 


Q 2 (1 + x 2 - 2£ 2 ) 
/i 2 (1 -x)(l - £ 2 ) 


0101 - 


(2.13) 


where Cf — (A 9 ,0 1 ^ for SU(N). We have cut off the transverse momentum integral at some 
scale Q and fi is the factorization scale separating hard and soft dynamics [18]. Also, we 
have taken | A A | to be small. For convenience, we take A 2 = 0. It is important to note that 
the entire a s dependency in Eq. (2.13) comes from the state and the operator is independent 
of interaction. The single particle matrix element receives a contribution upto 0(ck s ) from 

the normalization of the state. Taking into account the normalization contribution, we get 
i 


1 Here is the usual (principal value) plus prescription. 
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(2.14) 


F + 


1-C 2 5(1 



x) + 


(i+x 2 ~2a . 

(i-x) + (i-a ;j ‘ 


The end point singularity at x — 1 is canceled by the contribution from the normalization 
of the state to the single particle matrix element, similar to the hclicity independent case 
[23]. The splitting function can be easily extracted from the above expression : 


P qq = c, 


1 + x 2 - 2£ 2 

(i-x) + (i-e)' 


(2.15) 


This agrees with the known result of [26] (when replacing £ in Ref. [26] by 2£). 

Turning next to the hclicity flip part of the matrix element, we find that it solely arises 
from the mass term in the expression, Eq. (2.12), of the two-particle wave function . The 
form of the wave function shows that this contribution is suppressed. 

The hclicity dependent off-forward matrix element is conventionally parametrized in 
terms of the generalized quark distributions, 


fp+ _ ___ 


Uy(P') \H q (x,€,t) 7 + 7 5 + E q {x,i,ty ^- 


Th (Pi 


(2.16) 


where U\(P) is the quark spinor in our case. The light-front spinors for longitudinally 
polarized quarks are given in the appendix. Using Eq. (A3), and also the fact that the linear 
mass-dependent hclicity flip terms give suppressed contribution to the matrix element, we 
obtain that, E is suppressed (it has no logarithmic divergent part), provided A 2 - is small. 
We therefore get, 


H(x,U) = o 


,, lr f/1 . « s „ , Q 2 (3 , (l + x 2 -2£ 2 ) n 

i) = j [6(1 ~ *) + ^C,lo a -(-6( 1 - x) + (1 ^ )+(1 _ a ) 


The forward limit is easily obtained by putting £ = 0: 

1 oi x „ , Q 2 /3 


a 

2n 


(1 + x 2 ) 


H(x, 0,0) = - [5(1 -x) + ^ C f log ^ (-5(1 - x) + 


/x- 


(1 — x). 


(2.17) 


(2.18) 


The above expression can be identified with g\(x) for a dressed quark target, as calculated 
in [17]. This gives the intrinsic hclicity distribution for a quark dressed with a gluon in 
perturbation theory. 


B. Gluon Distribution 


In this section, we calculate the gluon distribution, 
p+ _ " / j~-AP+ 

where 

s +1 ~ 2 = 2. 


f dz-ei^PX | F + “(-y)U(y) I ™>- 


r ~ 2 t r P a) 


(2.19) 


( 2 . 20 ) 


6 



We use light-front gauge A + = 0. 

The Fock space expansion of the relevant part of the operator is given by, 

n. = ——— 

7r)3)2 

Here, A is the gluon helicity. We calculate the matrix element for a quark state dressed with 
a gluon. The Fock space expansion of the state is given by Eq. (2.7). The single particle 
sector does not contribute to the matrix element and the only contribution comes only from 
the two particle sector. 

The matrix element is given by, 

K = \ £ A / q± + (i_^ A± )^ 2 -^ 2 - ( 2 - 22 ) 

We have suppressed the quark helicity dependence of the wave functions and the sum over 
them. Using the full form of the two particle wave function, we find that the helicity flip 
terms proportional to the quark mass give suppressed contribution and the helicity non-flip 
part is given by, 



5] A J dkfd 2 ki J dk+d 2 kia\k i, A )a(k 2 , X)S(2xP + - kf - fc+). (2.21) 



x 2i r 


Cflog 



(1 — x) 2 - 


(2.23) 


where the first (second) term in the square bracket comes from the state with gluon helicity 
+ 1(—1). So we have, 




9 2n J * /i 2 x\J 1 — £ 2 


Using the parametrization of F + in terms of H„ and E„, one can write, 


1 _ r — r -A A 

EgX'x = (.P ') [H g (x, e, 07 + 7 5 + E g (x, C, t) 7 


2 M 


Ux(P). 


(2.24) 


(2.25) 


The fact that the helicity flip part of the matrix element is suppressed means that E g is also 
suppressed. So we get, 

Q 2 (i-(i-x) 2 -e) 


a, 


H g (x,£,t) = --Cflog 


2n 


k 


x{i-k, 2 ) 


(2.26) 


The splitting function can easily be extracted and is given by, 

99 f x(i - e) 

which again agrees with [26] when making the replacement of £ of [26] by 2£. Also, in the 
forward limit, Eq. (2.26) gives, 


(2.27) 


or, 


£r /- n m , Q 2 (!-(!- X) 2 ) 

H g (x,0,0) — 2^0 flog ^2 y 

(2.28) 

a O 2 

Hg{ 1 - X, 0, 0) = -kCflog-Xj(l + x). 

2n fi* 

(2.29) 


This gives the gluon intrinsic helicity distribution for a dressed quark target. In Eq. (2.29), 
we have taken 1 — x as the momentum fraction of the gluon, in order to compare with [17]. 
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III. TWIST-THREE DISTRIBUTION 


We now calculate the twist-three (transverse) component of the hclicity dependent off- 
forward distribution in perturbation theory. The matrix element of the transverse component 
is given by: 

Pyx = f ^ei^iP'X \ q-y | P\). (3.1) 

We calculate the above matrix element for a transversely polarized dressed quark state. 
As before, we work in light-front gauge A + = 0. The bilocal operator in this case can be 
written as, 

0 ±5 = $(-y )tW( : y) = y)aW”(y) + ^ _t (-y)«V^ + ( y)- (3-2) 

The operator involves the constrained field ^“(^-) and therefore it is called higher twist. In 
the light-front gauge, can be eliminated using the constraint equation, 

= -^-[a ± - {id 1 - + gA L ) + 7 °m]-0 + , (3.3) 

where the operator -qr is defined as, using antisymmetric boundary condition, 

y f(x~) = \f_ oo d y~< x ~ ~ y~)f(y~)- (3-4) 

The operator, in terms of the dynamical fields, can be written as, 

o ±5 = oi + ol + oy (3.5) 

where, 

= m4t i5T 4 + m (i^r 4t )‘ 1, < 3 - 6 ) 

c& = t'(-^)(-a 2 + + [ ( a 2 + )]* y. ( 3 . 7 ) 

o; = l 2 + a 3 A')<i(^) + g [-U**(-V)(-i^ 2 + <73^)]*( ( 3 - 8 ) 

Here, <f> is the two component fermion held, 

V"= g ■ (3-9) 

The Fock space expansion of $ is given by Eq. (2.3), with 77 being the two-component 
spinor. The operator has three parts : O ^ is the quark mass contribution, 0£ ± is the 

8 



quark transverse momentum contribution and Og is the quark-gluon interaction effect. The 
light-front expression clearly shows each contribution separately in light-front gauge. 

The longitudinally polarized dressed quark state is given in Eq. (2.7). The transversely 
polarized state is expressed in terms of the helicity states as, 

| k + , k\ s 1 } = -J=(| k + , k x , t)± | k + , k L , j», (3.10) 

with s 1 = ±m-R, where uir is the renormalized mass of the quark. Without any loss of 
generality, we take the state to be polarized along the x direction. 

The contributions to the matrix element coming from the three parts of the operator are 
given by, 


F 1 = 




x)^i 


x 


X 2 




,x-£ 


e 


,Q~ 


1 — x . 

+ 




(3.11) 


We have suppressed the quark helicity dependence of the wave function. Using the explicit 
form of the two particle wave function, 


Fl 


m 1 

p+ 


'ipl'ipi ^(i 


x) + 



2x(x — 2£ 2 ) y 
(l-x)(x 2 -^ 2 )^’ 


(3.12) 


= - i E / 


,x 




E AA 


i-A 

■ X-Z 

i-e 


1 -X , ,x + £ u q 2 £ 

+ - —- A -i<i 


,Q ± + 


-e 

1 — X 


i + C 


P+ X 2 - £ 2 


+ 


^)X ] J~F-Xa'^2\Y^ 




x 


— 9 

x z 


(3.13) 


This gives, 


Fl 


fcj- 


rn 1 r Q 2 a s (1 - x)(x 2 + £ 2 + 2x£ 2 ) 

p+^r^e G/ (x 2 -£ 2 )(i-£ 2 ) 


(3.14) 


The interaction part gives overlap contribution in terms of two-and one particle wave function 
and is given by, 


F 1 

9 


Cflog 


Q 2 a s m 1 
~^2^2P + 


5(1 — x) 


(3.15) 


As before, we have taken A 1 to be small. The interaction gives contribution only at the 
end point x — 1. Considering the normalization contribution to the single particle matrix 
element, we get, the total contribution, 


F 1 


m 1 


5(1 — x) + Cflog 


Q 2 a s 

H 2 2-7T 


(25(1 - x) + 


1 + 2a; (1 — £ 2 ) — x 2 y 

(i-x) + (i-a 


(3.16) 


Here, we have also considered the contribution of the normalization condition to the single 
particle matrix element, which cancels the end point singularity, similar to the twist-two 
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case. Here, m is the bare quark mass. The above expression has no singularity at x — £. 
In light-front theory, the linear mass term appearing in the light-front QCD Hamiltonian is 
renormalized as [27], 


rrin = m 


0 



(3.17) 


Here m R is the renormalized mass of the quark. The linear mass terms in light-front QCD 
Hamiltonian are associated with explicit chiral symmetry breaking [28]. Also, from the above 
results, we find that the three contributions, including the quark transverse momentum effect 
and the quark gluon interaction effect are proportional to the quark mass, which shows that 
the twist-three distribution is directly related to the dynamical effect of chiral symmetry 
breaking. In terms of the renormalized mass, we get, 


F 1 


m R 1 


5(1 — x) + Cflog 


Q 2 a s 
fi 2 2n 



1 + 2x(l — £ 2 ) — a: 2 y 

(l-u+a-a »’ 


(3.18) 


In the forward limit, this gives, 


F 1 


mR 

P+ 


5(1 — x) + Cflog 


Q 2 a s 

g 2 2n 



x) + 


(1 + 2x — X 2 )y 

(!-*)+ H ' 


(3.19) 


By comparing the rhs of the above equation with the transversely polarized structure func¬ 
tion gx for a dressed quark target [18], one obtains that, 


_ 2rriR _ 2S T 

F — ——gr — -rrrdTi 


P+ 


P+ 


(3.20) 


Since for a transversely polarized dressed quark, m R = St (see appendix). 


IV. EXAMINATION OF THE WANDZURA-WILCZEK RELATION IN 

PERTURBATION THEORY 


The twist-three matrix element is parametrized as [29], 


F ± = J-C 7 (P / ) ( 7 X 7 5 H + —7 5 E + — Y - Gj + jWG 2 + A ± I^7 5 G 3 
P+ v ' 2 M 2 M P+ 

+ ie ±v A v t-G,)u{P). 


(4.1) 


The light-front spinors for a transversely polarized quark are given by Eq. (A4) in the 
appendix. Using Eq. (4.1),(A5),(A6) we get, 

^ 1 = 7l + ( 4 - 2 ) 

which in the forward limit gives gr, since H(x, 0,0) = g\(x) and G 2 (x, 0,0) = g 2 {x) in 
the forward limit. Comparing with the result in the previous section, we see that Eq. (4.2) 
is in agreement with our result for a dressed quark. 
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Using Eq. (2.17) and (3.18) we get, 

h (- c l ni Q 2(X s\ 2(x + ^ 2 )n 

G 2 (x,£,t) = -C f log——[-5{l-x) + j, (4.3) 

which in the forward limit gives, 

~ 1 Q ^ CX r i 

G 2 (x,0,0) = -Cflog-^-^- -8(l-x)+2x. (4.4) 

2 g- 2 tt 1 

The above expression agrees with g 2 for a transversely polarized dressed quark target [18]. 

In the Wandzura-Wilczek approximation, where the quark mass as well as the quark- 
gluon interaction terms are neglected, the twist-three matrix element is given in terms of 
twist-two matrix elements as [12], 

= U(P')[^E(x,0 - ^ry+^H(x,0]U(P) 

+ j ^duG^{u,£)W + (x,u,£) + ie ±llk J ^duG k (u,£)W_(x,u,£), (4.5) 

where, 

a/ x 1/9 r) 

G“ («,?) = U(P')[tl(H + E)(u,Q + + tg^)E(u,0 

A /J. O O 

" 2F t+ ( u Tu + ? a?) + E) (u ’ 4 >1 U(p) ’ (4 ' 6) 


+ 4^(1 + u2- + (2P)E(u,o 


A/X CT CT 

- -^l + l8 u a7, + 


2£ ' ' ^ ft, ’ 

W±(x,u, £) are the Wandzura-Wilczek kernels given by, 


-.(x,u,€) = [e(x > O-r—7u - < O 7--77 

L («-0 (»-0 

,1/1/ . n 6>(u>x) . ..0(u 


6(u > x) .. 6{u < oh 

; -rr - 0 * < -0 7 , 

(m + 0 (« + 0 


(4.7) 


(4.8) 


Using the light-front spinors given in the appendix, we find, for a; > £ and = A 1 ; A 2 = 0 
that the WW relation for reduces to, 

2mR f , ^9(x £)0(u — x) . 

Ulw = J duH(u,0 y t (4.9) 

In the forward limit, the rhs becomes which gives the well known Wandzura- 

Wilczek relation for the transversely polarized DIS structure function g?. The twist three 
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vector distribution F 1 - is similarly expressed in terms of a Wandzura-Wilczek relation, how¬ 
ever it vanishes in the forward limit. Using the expression for H for a massive dressed quark 
in perturbation theory, 


fp i 

r ww 


pvi-d (l " 0 qiw + k 
+ (T^a ( (1 +o^Arf) - 1 +d 


(4.10) 


Comparing the above expression with Eq. (3.18), we see that the WW relation is not satisfied 
for a dressed quark state in perturbation theory, similar to the forward case. This is not 
surprising because in the WW approximation, the mass of the quark as well as the explicit 
interaction dependence of the operator is neglected, whereas, we have obtained the full result 
in perturbative QCD for a massive quark.The effect of quark transverse momentum in gx 
was investigated in a covariant parton model approach in [30]. Also, it is known that in the 
forward limit, the WW relation is violated in perturbation theory [18]. However, BC sum 
rule is satisfied [31,18,32], 

The quark mass effect can be incorporated in the derivation of the off-forward WW 
relation [29]. This gives an additional contribution to F L which is of the form (see Appendix 

B), 


mass j^+ 


for £ < x < 1, where 


x 


x 2 — £ 2 


/ x (x,£,A)+/ dy 


y 2 +e 
* ~“{y 2 -e ) 2 


/ ± (2/^, A ) 


/ ± (z,£, A) = 




,Z 


2J 2vr 

We use the parametrization [34], 


e 3 ‘ ' ‘(f'a' I V>(-y)'» +± 7 5 </’(y) I PA >- 


(4.11) 


(4.12) 


\ I | &-y)«7 + W(y) I p A> 


1 _ r ~ 7 /\ P/j /\ ry „ 

= A') + HI M2 a 13 + Eh 2M ^ 


+ E q r 


T M 2 


M 


U(P, A). 


(4.13) 


Here M is the mass of the state. We calculate the above matrix element for a transversely 
polarized dressed quark state in perturbation theory. Using the relations of light cone 
spinors, and also using the normalization of the transversely polarized state, we obtain, 


Ht — ~ 1^(1 ~ x) + CflogQ-—(-S(l — x) + 

T 2^ v ; f y /i 2 27r A 2 v ; (l-x)+(l-£ 2 ) 


2(x-a 


which in the forward limit gives h\{x) for a dressed quark : 

'3 


2x 


hl (x) = i[5(l - *) + C,tog^ty( 1 - *) + (1 _ 


(4.14) 


(4.15) 
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Next, we investigate the mass corrections to the WW relation and the ‘genuine twist three 
contribution ’ to the matrix element in somewhat more detail. In the forward limit, F 1 
corresponds to Qt- We can write, in the forward limit, 


/ n f 1 , 9i(y) , m Mx) r 1 fti(2/)\ 

9 t\X) = / dy - + — ( - / dy — —)+g y T (x) 

Jx y M v x Jx y 2 ' 


(4.16) 


where g^(x) is the so-called ‘genuine twist three’ contribution to gr- If we neglect this, we 
get the WW relation with the quark mass correction, 


/ N f 1 r 9i(y) , m (hi{x) 

9t(x) = L dy ^r + m(— 



(4.17) 


Here m is the quark mass and M is the mass of the target. It is very important to note that 
in perturbative calculation, m has to be renormalized. Taking n-tli moment of both sides 
of Eq. (4.16) we get, 


92 = 


n 


+ 1 


9 1 + 


M n + 1 


«i 


T HT, 


(4.18) 


where a n = /g 1 dxx n a(x). Using the expressions of gi(x),g 2 (x) and h\{x), the moments can 
be directly calculated, 


1 a , 


92 2 2 ? y 2 


Cflog 


Q 2 


n 


n + T 


(4.19) 




1 V 1 )! 

(n + l)(n + 2) 2 pzy- 


(4.20) 


h 


n 

1 


1 r 1 o: s 

1 + 2i C/l ° 9 ^ 


2 L 


Q 2 f 3 



(4.21) 


Using these and also the renormalization of quark mass given by Eq. (3.17), we obtain, 


9t q 


la* 

2 2vr 


Cflog 



n n (3 

n+2^n+lW 


2n + 3 y 

(n + 1 )(n + 2 y-' 


(4.22) 


For n = 0, the rhs of the above equation gives zero, which proves the BC sum rule. For 
n = 1 , the rhs of Eq. (4.22) also yields zero, which gives the Efremov-Leader-Teryaev sum 
rule with the correction due to quark mass. 

Next, we use Eq. (4.16) to extract the ‘genuine twist three’ part of gr- 
The 0(a s ) part of Qt can be separated into two parts, 



9ta 


9tb- 


(4.23) 


Here g^\ is the WW part with the mass corrections and g^ B is the ‘genuine twist three 
part’. Using Eq. (4.16) and also the expressions of gi(x),hi(x) we get, 
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(4.24) 


9tb = 


1 Os 
22n 


Cfl ° 9 ~I [^( 1 logx 


It is interesting to compare the rhs of the above equation with the forward limit of Eq. 
(3.15). This shows that Eq.(3.15) does not give the full ‘genuine’ twist three contribution 
but only a part of it. Also from the above expression it is easy to check that the first and 
second moments of are zero. 

We stress that, the quark mass plays a very important role in the twist-three matrix 
element, and also, in our case, it is essential to obtain a transversely polarized state, since 
St = '/rift, the renormalized mass of the quark. Our result shows that the twist-three 
generalized distribution is directly related to the chiral symmetry breaking dynamics in 
light-front QCD. 


V. SUMMARY AND DISCUSSIONS 

To summarize, in this work, we have investigated the off-forward matrix elements of 
the light-cone axial vector operator. We have calculated the matrix elements of the plus 
and transverse components of the operator for a dressed quark in light-front Hamiltonian 
perturbation theory. This approach allows us to express the distributions in terms of light- 
front wave functions. We have restricted ourselves to the kinematical region £ < x < 1. 
In this case, the overlaps of three-particle and one-particle wave functions are absent. We 
obtained the splitting functions for the evolution of the helicity dependent twist two quark 
and gluon distributions in a straightforward way. We showed that the singularity at x — 1 
is canceled by the contribution from the normalization of the state, similar to the helicity 
independent case calculated earlier. The twist-two distributions reduce to the quark and 
gluon intrinsic helicity distributions for a dressed quark target in the forward limit. The 
twist-three distribution is expressed entirely in terms of the dynamical fields in the light- 
front gauge. This calculation shows that for the twist-two distributions, the entire interaction 
dependence comes from the state whereas the operator has free field structure, but in the 
case of twist-three, both the operator and the state introduce interaction dependence. The 
operator has three parts, explicit mass dependent term, quark-gluon interaction term and 
a term containing the quark transverse momentum effect. The calculation of this matrix 
element for a transversely polarized dressed quark shows that all the three contributions 
are proportional to the quark mass. Using the renormalized quark mass niR in light-front 
Hamiltonian perturbation theory, we found that in the forward limit, F x is proportional to 
St9t, where St is the transverse polarization of the state, St = rriR in our case. It is known 
that in light-front Hamiltonian QCD, chirality is the same as helicity, and the terms that 
cause helicity flip in the light-front QCD Hamiltonian are explicit chiral symmetry breaking 
terms. These terms are linear in the quark mass. It is interesting to note that the quadratic 
mass terms do not flip helicity, however, they are suppressed here. Therefore, we concluded 
that F 1 - is directly related to the chiral symmetry breaking dynamics in light-front QCD 
and the quark mass plays an important role. In particular, a finite mass is necessary to 
have a transversely polarized quark state. We have calculated the same off-forward matrix 
element in the Wandzura-Wilczek approximation and found that the actual result for a 
massive dressed quark deviates from the WW approximated form. The violation of the 
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WW relation for gx for a massive quark is known in perturbation theory and our result 
reduces to gx for a massive dressed quark in the forward limit. It is to be noted that 
in the case of nucleons, the quark intrinsic transverse momentum effects and the quark- 
gluon coupling dynamics play a more complicated role and the pure quark mass effects in 
perturbative QCD may be suppressed by where M is the hadron mass. We have also 
calculated the quark mass correction to the off forward WW relation which in the forward 
limit reduces to a term proportional to hi(x). We extracted the ’genuine twist three part’ 
of gx and showed that both BC and ELT sum rules are satisfied. 

It is known that in the kinematical region 0 < x < £, contribution comes from the overlap 
of three-and-one-particle wave functions. The GPDs in this region have a different type of 
evolution (Brodsky-Lepage). It will be interesting to investigate the GPDs for a dressed 
quark in this kinematical region using this approach and to check the various moment 
relations in the whole range of x, 0 < x < 1. Another interesting topic for future work is to 
investigate the A -1 dependence of the GPD’s in the frame £ = 0. 
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APPENDIX A: LIGHT-FRONT SPINORS 


The light-front spinors for longitudinally polarized quark of mass M and momentum P 
and hclicity up and down, respectively, are given by [35], 


U \{P) = 


1 


V2P+ 

Using these, we get, 


( P+ + M\ 
P 1 + iP 2 

1 

/ -P 1 + iP 2 \ 
P + + M 

P+ - M 

’ ^ ) V / 2 P+ 

P 1 - iP 2 

V P 1 + iP 2 j 


\ -P+ + M ) 


u^p'^^p) = 2^i~ ep 

2 iM 


>+ 


f/ T (P'bRR) = 




Also, 


£/ t (F') 7 + 7 5 C/i(f) = 0 

-A 1 +iA 2 


U^P')^U { (P) = 


VR 2 


(Al) 


(A2) 


(A3) 
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APPENDIX B: MASS TERM IN WW RELATION 

In this appendix, we give an outline of the derivation of Eq. (4.11). Using the approach 
described in [29], and taking into account the quark mass m, one gets, 

= j c ^-e-^(P'S' | q-i-)f>rV(y) I PS) 

= + (Bl) 

were A = \P + z~, M a (x, £,t) is the mass term and X a (x,^,t) are all the other terms 
considered in the WW approximation. Here we concentrate on the mass term given by [33], 

/ d\ r 1 

—e~ iXx \ / duu e i(1 -“ )?A + e - i(1 -“ )?A 
2-7T JO L J 

\(P'S' I | PS). (B2) 

If we define, 

/“(*,£, A) = 1 / I I PS), (B3) 

we get, from Eq. (B2), 
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777 f ) /*1 

M “ (x ^ t) = p^TxL du 


/“( 


x + (1 -u)f 


u 


,c,A) + r( 


® - (1 - u)f 


u 


,e,A) 


(B4) 


Changing the variable, y = ,l±(1 — we obtain, for x > £, 


M a (x, £, t) = 


2m 

P+ 


x 


x 2 -e ; 


-/ Q (x,£,A)+/ % 


7 2 + e 2 


- "(?/ 2 -a : 


A) 


(B5) 


17 



REFERENCES 


[1] D. Muller, D. Robaschik, B. Geyer, F. M. Dittes, J. Horejsi, Fortsch. Phys. 42, 101 
(1994). 

[2] X. Ji, J. Phys. G 24, 1181 (1998). 

[3] A. V. Radyushkin, hep-ph/0101225, published in ”At the Frontier of Particle 
Physics/Handbook of QCD”, ed. M. Shifman (World Scientific, Singapore, 2001). 

[4] K. Goeke, M. V. Polyakov, M. Vanderhaeghen, Prog. Part. Nncl. Phys. 47, 401 (2001). 

[5] S. J. Brodsky, M. Diehl, D. S. Hwang, Nncl. Phys. B 596, 99 (2001). 

[6] M. Diehl, T. Feldmann, R. Jakob, P. Kroll, Nncl. Phys. B 596, 33 (2001). 

[7] M. Diehl, T. Feldmann, R. Jakob, P. Kroll, Enr. Phys. J C8, 409 (1999); B. C. Tibnrzi 
and G. A. Miller, Phys. Rev. C 64, 065204 (2001); A. Mukherjee, I. V. Musatov, H. C. 
Pauli, A. V. Radyushkin, hep-ph/0205315. 

[8] A. Airapetian et al, (HERMES Collaboration), Phys. Rev. Lett. 87, 182001 (2001). 

[9] S. Stepanyan et al, (CLAS Collaboration), Phys. Rev. Lett. 87, 182002 (2001). 

[10] C. Adloff et al, (HI Collaboration), Phys. Lett. B 517, 47 (2001). 

[11] P. R. Saull et al, (ZELTS Collaboration), hep-ex/0003030. 

[12] A. V. Belitsky and D. Muller, Nucl. Phys. B589, 611 (2000); N. Kivel, M. V. Polyakov, 
A. Schafer and O. V. Teryaev, Phys. Lett. B497,73 (2001). 

[13] S. Wandzura and F. Wilczek, Phys. Lett. B72, 195 (1977). 

[14] P. L. Anthony et al, (E155 Collaboration), hep-ex/0204028. 

[15] W. M. Zhang and A. Harindranath, Phys. Rev. D48, 4881 (1993). 

[16] A. Harindranath, R. Kundu, W. M. Zhang, Phys. Rev D 59, 094013 (1999). 

[17] A. Harindranath and R. Kundu, Phys. Rev. D 59, 116013 (1999). 

[18] A. Harindranath and W. M. Zhang, Phys. Lett. B408, 347 (1997). 

[19] A. Harindranath, A. Mukherjee, R. Ratabole, Phys. Rev. D 63, 045006 (2001). 

[20] A. Harindranath, R. Kundu, A. Mukherjee and J. P. Vary, Phys. Lett. B 417, 361 
(1997); Phys. Rev. D 58, 114022 (1998). 

[21] A. Mukherjee and D. Chakrabarti, Phys. Lett. B 506, 283 (2001). 

[22] R. Kundu and A. Metz, Phys. Rev. D 65, 014009 (2002). 

[23] A. Mukherjee and M. Vanderhaeghen, Phys. Lett. B 542, 245 (2002). 

[24] H. Burkhardt and W. N. Cottingham, Ann. of Phys, 56, 453 (1970). 

[25] A. V. Efremov, O. V. Teryaev and E. Leader, Phys. Rev. D 55, 4307 (1997). 

[26] X. Ji, Phys. Rev. D55, 7114 (1997). 

[27] R. J. Perry, Phys. Lett. B 300, 8 (1993); A. Harindranath and W. M. Zhang, Phys. 
Rev. D 48, 4903 (1993). 

[28] K. G. Wilson, T. S. Walhout, A. Harindranath, W. M. Zhang, R. J. Perry and S. D. 
Glazek, Phys. Rev. D 49, 6720 (1994). 

[29] D. V. Kiptily and M. V. Polyakov, Preprint, hep-ph/0212372. 

[30] J. D. Jackson, G. G. Gross and R. G. Roberts, Phys. Lett. B 226, 159 (1989). 

[31] G. Altarelli, B. Lampe, P. Nason and G. Ridolh, Phys. Lett. B 334, 187 (1994). 

[32] M. Burkardt and Y. Koike, Nucl. Phys. B 632, 311 (2002). 

[33] D. V. Kiptily and M. V. Polyakov, private communication. 

[34] M. Diehl, Eur. Phys. J. C 19, 485 (2001). 

[35] A. Harindranath, An Introduction to the Light Front Dynamics for Pedestrians in Light- 


18 



front Quantization and Non-perturbative QCD , Ed. J. P. Vary and F. Wolz, published 
by Internatinal Institute of Theoretical and Applied Physics, Ames, Iowa, USA (1997). 
[36] A. Mukherjee, Ph. D. Thesis, University of Calcutta, hep-ph/0106167. 


19 



